Can the WMAP Haze really be a signature of annihilating neutralino dark matter? 
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Observations by the Wilkinson Microwave Anisotropy Probe (WMAP) satellite have identified 
an excess of microwave emission from the centre of the Milky Way. It has been suggested that 
this WMAP haze emission could potentially be synchrotron emission from relativistic electrons 
and positrons produced in the annihilations of one (or more) species of dark matter particles. 
In this paper we re-calculate the intensity and morphology of the WMAP haze using a multi- 
linear regression involving full-sky templates of the dominant forms of galactic foreground emission, 
using two different CMB sky signal estimators. The first estimator is a posterior mean CMB map, 
marginalized over a general foreground model using a Gibbs sampling technique, and the other is 
the ILC map produced by the WMAP team. Earlier analyses of the WMAP haze used the ILC 
map, which is more contaminated by galactic foregrounds than the Gibbs map. In either case, we 
re-confirm earlier results that a statistically significant residual emission remains after foreground 
subtraction that is concentrated around the galactic centre. However, we find that the significance of 
this emission can be significantly reduced by allowing for a subtle spatial variation in the frequency 
dependence of soft synchrotron emission in the inner and outer parts of the galaxy. We also re- 
investigate the prospect of a neutralino dark matter interpretation of the origin of the haze, and find 
that significant boosting in the dark matter annihilation rate is required, relative to that obtained 
with a smooth galactic dark matter distribution, in order to reproduce the inferred residual emission, 
contrary to that deduced in several recent studies. 

PACS numbers: 93.35. +d, 98.58.Ca, 98.80.-k, 41.60.Ap, 03.50.-z, 41.60.-m, 78.70.Ck, 98.70.Vc, 98.35.Jk 



I. INTRODUCTION 

As the statistical error bars on Cosmic Microwave 
Background (CMB) measurements shrink, the impor- 
tance of correct foreground subtraction grows. The sep- 
aration of the various forms of Galactic foreground emis- 
sion from the primordial CMB signal is both critical to 
correct estimation of the CMB statistics and very diffi- 
cult. The application of foreground subtraction to the 
all-sky microwave measurements from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) has been extensively 
studied jj]. 

In WMAP's frequency range 2Z < v < 94 GHz 
the dominant sources of foreground signal come from 
emission in the interstellar medium: free-free (thermal 
bremsstrahlung), dust and synchrotron emission. A cor- 
rect separation of these sources, which are the dominant 
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signal over large parts of the sky, allows both sensible use 
of the CMB data for cosmological parameter estimation 
and for an improved understanding of the emission mech- 
anisms in the ISM. The use of prior information about 
the likely spatial and frequency dependence of the sources 
makes this separation possible with current data. 

During the study of Galactic foreground emission 
evidence has emerged for an additional, possibly ex- 
otic, emission mechanism responsible for an excess of 
microwave radiation emerging from the Galactic cen- 
tre (GC) and distributed approximately symmetrically 
around it |2j. This currently unexplained excess emis- 
sion has been referred to as the WMAP haze, and its 
origin has been the subject of much debate. 

The WMAP Haze was originally suggested to be an 
excess in free-free emission, originating from hot ionised 
gas with temperature T ^ 10^ K, owing to its al- 
leged hard spectrum [3 . However, such gas is ther- 
mally unstable, and with an insufficient abundance of 
gas at hotter or cooler temperatures, evident from con- 
straints on the intensity of recombination lines/x-rays, 
such an explanation of the signal seems unlikely ^ |4]. 
The hard spectrum has also led to speculation (see 
e.g. [SI El |7l E]) that the origin of the haze is more 
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likely to be hard synchrotron foreground emission result- 
ing from a source distinct from that responsible for the 
softer synchrotron emission commonly associated with 
Galactic foregrounds. An astrophysical origin of such a 
hard synchrotron spectrum, such as that produced from a 
distribution of ultra-relativistic electrons, possibly result- 
ing from supernovae explosions, with a flux large enough 
to account for the haze, seems difficult to realise [8 . 

Here, we re-evaluate the recent proposal that the 
WMAP haze may originate from synchrotron emission 
produced by a distribution of high-energy positrons and 
electrons resulting from the self-annihilation of super- 
symmetric neutralino dark matter. What differentiates 
our study from those past [2 O El [3 is that in addition 
to using a standard set of templates, designed to model 
Galactic foreground emission, in a template-based, multi- 
linear fit to the CMB-subtracted WMAP data, we use a 
by-product of the Gibbs sampling method of CMB pa- 
rameter estimation as our CMB contribution to the data 
[10], and compare our results with those calculated us- 
ing the internal linear combination method (ILC) utilised 
in previous work [2^. Gibbs sampling is a Monte-Carlo 
Markov Chain method that generates posterior samples 
of the signal map, power spectrum, and foreground com- 
ponents. The map from this method is the output of a 
well-understood statistical process that is known to pro- 
duce good results, and may be one of the cleanest CMB 
temperature maps currently available. 

There have been many other proposals for the origin of 
the WMAP haze involving dark matter candidates other 
than supersymmetric neutralino dark matter. These in- 
clude, exciting dark matter (XDM) particles [7 which an- 
nihilate to produce a hard spectrum of centred around 
the GC. Also, compact composite object (CCO) dark 
matter [TT], which are macroscopic assemblies of quarks 
in a colour superconducting state whose finite electric 
charge is neutralised by an atmosphere of bound to 
each CCO. It is then proposed that these annihilate 
with matter in the ISM to produce radiation which con- 
tributes to a number of unexplained astrophysical ob- 
servations, including the WMAP Haze and the 511 keV 
excess from the GC (however the prospects for a CCO 
explanation of the 511 keV signal are considered by some 
to be unlikely ^2 )• Also, Sommerfeld enhanced dark 
matter has been proposed to explain the WMAP haze 
[13 , where the annihilation rate of TeV-mass particles 
are significantly boosted in low velocity environments, 
owing to an additional Yukawa-like attractive force be- 
tween colliding particles. 



II. WMAP DATA AND TEMPLATES FOR 
GALACTIC FOREGROUNDS 

The three forms of Galactic foreground emission that 
dominate in the range of frequencies observed by WMAP 
are: free-free (or thermal bremsstrahlung) emission from 
hot ionised gas; dust emission, in two forms: firstly 



from thermal emissions from dust particles interacting 
with the ambient radiation field, and secondly from spin- 
ning dust particles excited by collisions with ions travers- 
ing the ISM [TT, and synchrotron emission originat- 
ing from electrons which have been shock-accelerated to 
GeV-energies by supernova explosions. Since microwave 
radiation is not significantly attenuated by the parti- 
cles involved in these emission mechanisms, the all-sky 
maps produced from the CMB-subtracted WMAP data 
should theoretically be reproducible from a linear com- 
bination of templates modelling the various foreground 
components, constructed from independent data maps 
of the three dominant components of foreground emis- 
sion. But firstly, we briefly discuss the nature of the 
CMB-subtracted WMAP data that is to be used in our 
anlaysis. 



A. CMB-subtracted WMAP data 

We utilise the Syr WMAP data maps [T5^ provided 
in the HEALPix [16^ pixelization scheme. We construct 
temperature maps at a given frequency using a weighted 
average of all the WMAP differencing array (DA) maps 
at that frequency. The weighting, which is different at 
each pixel p, is given by 1/cr^^, where dpi = (Ji/^/n^ is 
the error in the pixel for the given DA, where is the 
characteristic noise of the DA and Up is the number of 
times the pixel is measured. Such data has been released 
by WMAP [17 . Our CMB map ^ is the output of 
the Gibbs sampling method described, implemented and 
applied to the 3-year WMAP data in [TO . This CMB 
map was produced by marginalizing the joint foreground- 
CMB posterior over a free low-frequency foreground am- 
plitude and spectral index at each pixel, as well as a 
single amplitude relative to the dust template described 
below [19 . These additional degrees of freedom allows 
for a much more accurate CMB reconstruction than that 
provided by the WMAP ILC map. However, all results 
presented in this paper are computed for both the Gibbs 
and the ILC maps, to assess the impact of residual fore- 
grounds on our results. In general, we find only small 
differences between the two data sets, and no main con- 
clusions depend significantly on the CMB map choice. 

We used the point source mask recently updated and 
distributed with the latest 5-year data release by WMAP 
[20 to mask areas of the sky which may contaminate our 
fit.^ 



B. Free- free template 

Free-free (or thermal bremsstrahlung) emission origi- 
nates from coulomb interactions between free electrons 
and hot interstellar gas. Since free-free emission is pro- 
portional to the emission measure EM= J nidi, where 
He is the spatially-dependent number density of inter- 
stellar electrons, along a given line of sight, maps of Ha 
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recombination line emission (which is also proportional 
to the emission measure) can be used to approximately 
trace the morphology of free-free emission in the Galaxy. 
In fact, Ha maps at intermediate and high latitudes are 
the only effective free- free templates available. In this 
study we use the all-sky Ra template described in [21] 
which is a composite of data obtained from the WHAM 
Fabry-Perot survey of the northern sky [22 (which pro- 
vides a good distinction between the Galactic emission 
and the geocoronal Ha emission, which is caused by so- 
lar Ly-f] excitation of HI in the exosphere) the SHASSA 
filter survey of the southern sky |23j and the variable 
resolution VTSS filter survey of the northern hemisphere 
[24] . When using the Ra map as a template for free- free 
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FIG. 1: The Free-free template map, masked using the 
WMAP five-year point source mask described in § |II A| and 
presented in units of log[Rayleighs(i?)], where IR = 10^/47r 
photons cm~^ s~^ sr~^. 

emission, it is necessary to correct for extinction by dust 
(both in front of and mixed with the warm gas gener- 
ating the Ha fiux). We correct our template for dust 
extinction by following the treatment in ^ [21j : we as- 
sume that the warm gas is uniformly mixed with dust in 
a cloud possessing optical depth r, the observed spectral 
intensity Iu,ohs. is then related to the emitted intensity 

L,em. by 

/.,obs. = ^(l-e--), (1) 

where r is computed by multiplying the all-sky 100 /im 
dust map (here, in units of mK antenna temperature) 
produced by Schlegel, Finkbeiner and Davis (SFD hence- 
forth) [25^, in units of E{B — V) magnitudes reddening, 
by a factor 2.65/1.086 [21]. We then further reduce the 
emitted fiux by masking out all regions where the dust 
extinction A{Ra) = 2mE{B - V) > 1. 

The conversion of dust-corrected Ha intensities to EM 
and then to free-free emission intensities is well under- 
stood. The corresponding free-free brightness tempera- 
ture T5 is then related to the electron temperature and 
EM by 

n cx T-^-^^y-'^-^^ X EM, (2) 

where over the range of frequencies v observed by 
WMAP, Te is approximately constant (c^ 8000K)[4l[26]. 



The resulting template map is displayed in fig. [T] 

C. Dust template 

Thermal dust emission is produced by microscopic in- 
terstellar dust grains vibrating in thermal equilibrium 
with the ambient radiation field. Here, we use the tem- 
plate map of dust emission by Finkbeiner, Davis and 
Schlegel (FDS henceforth), evaluated at 94 GHz [19 . It 
has also been proposed in [^ that such a map should 
closely trace the electric dipole emission expected to be 
produced by the smallest dust grains, which have a finite 
dipole moment, excited into rotational modes mainly by 
collisions with surrounding ions [28] [22 . The dust tem- 




FIG. 2: The dust template map, masked using the WMAP 
five-year point source mask, and presented using arbitrary 
units. 

plate map is displayed in fig. [2j 

D. Synchrotron template 

Galactic synchrotron emission is mainly produced by 
electrons in the vicinity of supernovae explosions which 
are shock-accelerated to relativistic (at least MeV) ener- 
gies and then subsequently lose energy, mainly as inverse 
Compton radiation from scattering with CMB photons 
and starlight, and synchrotron radiation from their in- 
teraction with the Galactic magnetic field. This emis- 
sion is best measured at low frequencies (i.e. < IGHz ) 
where contamination from other sources, mainly free- free 
emission, is its lowest. A full-sky map of this emission 
was measured at 408 MHz by Haslam et al. [30 , which 
we use here as our template of Galactic synchrotron fore- 
ground emission. This synchrotron template is displayed 
in fig. [3| 

The spectral dependence of the synchrotron brightness 
temperature, T5 oc z/^, is expected to be highly spa- 
tially dependent, especially (hardening) near to recent 
supernova activity. Past analyses of low frequency maps 
of Galactic synchrotron radiation have deduced wide- 
ranging values of the index jS consistent with the available 
data, with a value in the vicinity of -3 [H |3Tj . 



FIG. 3: The synchrotron template map, masked using the 
WMAP five-year point source mask, and presented in units 
\og[T{K)], where T is antenna temperature. 



E. Template fitting procedure 

We perform a fit of the CMB-subtracted WMAP data 
described in § |II A| using maps in all five WMAP fre- 
quency bands (K, Ka, Q, V and W, corresponding to 
frequencies of 22.8, 33.0, 40.7, 60.8 and 93.5 GHZ respec- 
tively) by using a template-based multi-linear regression. 
We follow a treatment identical to that described in |9] 
and find the solution vector a for the template coefficients 
that minimizes the of the linear equation 



Pa = w, 



(3) 



where P describes the templates, w the data and a the 
coefficients of the templates in each frequency band. 

We work with maps all downgraded to resolution 
^side = 64 and smoothed using a beam of FWHM=3°. 
We extend the WMAP 5 year mask to cover the region 
where edge effects in the smoothing come into play. We 
also manually remove point sources that appear as strong 
residuals during the analysis. 



III. 3-TEMPLATE FIT 

In this section we demonstrate our template-based 
multi-linear fitting procedure by performing a basic linear 
regression involving the full-sky free- free, dust and syn- 
chrotron templates described and displayed in the previ- 
ous section. 

We characterise the quality of our fits by using the 
statistic Xred ~ where is the Euclidean norm of 

the residual vector r = Pa — w, divided by the WMAP 
la (noise) error maps a 
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and u = Np — Na^ where Np is the number of unmasked 
pixels in each map and Na is the number of elements of 
a. 





FIG. 4: Residual emission, in units of mK antenna tempera- 
ture, in the K (top), Ka (middle) and Q (bottom) frequency 
bands, when using the Gibbs (left) and ILC (right) CMB es- 
timators. 



In fig.[4]we see the result of this procedure, which shows 
the residual emission in the K (top), Ka (middle) and Q 
(bottom) WMAP bands when using the Gibbs (left) and 
ILC (right) CMB estimators. (We do not display residual 
maps for the V and W bands since the large experimental 
errors aij heavily suppress the variance at these frequen- 
cies.) The overall Xred 

of the fit is equal to 7.85 (11.05) 
for the Gibbs (ILC) analysis. To further gauge the signif- 
icance of the residual emission we calculate the following 
parameter 



1 



1 



(5) 



where Xij and a, 



i^j are the respective values of the jth 
pixel in the zth frequency band associated with the resid- 
ual map and CMB la error map associated with that 
band. {x)ij = is the expectation value of x^^j, since 
we expect, as a first approximation, that we will obtain 
a perfect fit to the data. Therefore, looking at the form 
of ([5| we can say that regions of statistically significant 
residual signal possess values C.'f^j ^ 1. 

In fig. [5] we display maps of in the K (top), Ka 
(middle) and Q (bottom) bands, for both Gibbs (left) 
and ILC (right) analyses, truncating the maps with a 
minimum of Cf ^=l, i-e. only displaying regions of sta- 
tistically significant residual emission. The values of Q 
corresponding to the residual maps in each of the five 
frequency bands are: 5.54 (6.59), 0.88 (1.45), 1.08 (2.12), 
0.089 (0.174) and 0.028 (0.271) for the K, Ka, Q, V and 
W bands respectively, for the Gibbs (ILC) analysis. The 
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FIG. 5: Maps of the parameter (^fj in the K (top), Ka (mid- 
dle) and Q (bottom) frequency bands, when using the Gibbs 
(left) and ILC (right) CMB estimators, where we truncate the 
maps using a minimum of (fj = 1. 
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FIG. 6: Spectral dependences of best fitted free- free (blue 
markers), dust (red markers) and synchrotron emission (green 
markers) foreground components resulting from our fit when 
using the Gibbs (circles) and ILC (crosses) CMB estimators. 
Also plotted are several lines of best fit associated with the 
best fit spectral indices for the data points (dashed black 
(Gibbs) and dashed magenta (ILC) ) associated with free- 
free and synchrotron components, and also for dust emission 
in the K-Ka and V-W dust frequency bands. For clarity, we 
have slightly displaced the frequency of the two sets of data 
markers associated the Gibbs and ILC CMB estimators by 
Alog]^Q(z//GHz)=+0.5 and —0.5 respectively. 



1(7 errors associated with these values are all of order 
10~^. Hence it is clear that a statistically significant (i.e. 

> 1) residual emission remains, at least for the K 
band, which appears to be concentrated mostly around 
the GC. In fact, the respective values of calculated 
solely for the sky region within 50° of the GC increase 
to 14.69 (16.59), 1.65 (2.42), 1.60 (2.84), 0.11 (0.20) and 
0.027 (0.30). The la errors associated with these values 
are all of the order 10~^ or less. The errors associated 
with the results presented throughout this paper are la 
errors (unless stated otherwise). They were calculated 
using a random sample of 1000 realizations of the CMB 
data maps, constructed using the average CMB-removed 
WMAP data and the la error maps a. 

Using the elements of the solution vector a we can cal- 
culate the spectral dependence of the various foreground 
components. In fig. [6] we plot the spectral dependences 
of the free- free (blue markers), dust (red markers) and 
synchrotron emission (green markers) resulting from our 
fit using the Gibbs (circles) and ILC (crosses) CMB es- 
timators. We obtain best fit values of spectral index 
a equal to -1.93 (-1.80) and -2.93 (-3.13) for the free- 
free and synchrotron emission components respectively, 
when using the Gibbs (ILC) estimators, where we as- 
sume that the foreground intensity 1^ scales with fre- 
quency as (X These results are consistent with 
those determined in previous studies discussed in § II 
and §|IID[ To further illustrate this we have plotted lines 



of best fit with these gradients (dashed black (Gibbs) and 
dashed magenta (ILC) ). For clarity, we have slightly 
displaced the frequency of the two sets of data mark- 



ers associated the Gibbs and ILC CMB estimators by 
Alog;Lo(^/GIIz)=+0.5 and —0.5 respectively. 

Unlike free-free and synchrotron emission, the emission 
from dust does not have a constant spectral index, but 
varies from -2.93 (-2.99) between the K and Ka bands 
to +1.14 (+1.23) between the V and W bands when us- 
ing the Gibbs (ILC) analysis, with a minimum intensity 
occurring around 50 GHz. Again we illustrate this by 
plotting lines of best-fit with these gradients. Such re- 
sults are consistent with the notion that there are two 
components to dust emission, dominating at the opposite 
extremities of the WMAP frequency range. It has been 
suggested that emission from both spinning dust, domi- 
nating in the K, Ka and Q bands with a large negative 
index of approximately -2.85, and thermal dust, domi- 
nating in the W band with a positive spectral index of 
approximately +1.7, may be consistent with the observed 
spectral behaviour of the dust emission (see e.g. [27|l32]), 
however we do not to pursue the investigation of these 
topics further. 

We now briefly consider the correlation between the 
residual haze and the various foreground components. 
We define the elements of correlation matrix, 0x,Y7 in- 
volving foreground template X and residual or fore- 
ground template Y (with elements Xi and yi respectively) 



as 
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TABLE I: Correlation matrix involving the free-free, dust and 
synchrotron foreground templates together with the map of 
residual emission in the K-band following the 3-template fit. 



In table [l| we display the correlation matrix involving each 
of the three foreground templates and the K-band resid- 
ual map when using the Gibbs CMB estimator (the corre- 
lation matrix when using the ILC estimator is extremely 
similar). One can clearly observe that the haze emission 
in the K-band is significantly more correlated with syn- 
chrotron emission (0.48) than either the dust (0.19) or 
free- free (0.11) emission. We illustrate this correlation in 
the following section where we demonstrate how we can 
substantially improve the quality of the fit to the data by 
allowing the spectral dependence of the soft synchrotron 
emission for regions close to the GC to vary indepen- 
dently from that associated with regions far away from 
the GC. 



IV. FIT USING MULTIPLE SYNCHROTRON 
TEMPLATES 

In this section we investigate whether the significant 
correlation between the residual haze and the soft syn- 
chrotron galactic foreground emission demonstrated in 
section §[111] can be explained by a subtle spatial varia- 
tion in the spectral dependence of soft synchrotron galac- 
tic foreground emission. Here, we adopt a simple, yet 
illustrative model to allow for such spatial variation by 
splitting the synchrotron template into two regions (one 
central and one peripheral to the GC), and allow the 
two to have different spectral indices. We propose, for 
simplicity, that the boundary separating the two regions 
should be symmetrical around the GC, corresponding to 
all lines of sight with inclination relative to the GC. 

In fig. [7| we illustrate an example case where the 
Haslam et al. synchrotron template is segregated into 
central and peripheral regions according to the prescrip- 
tion described above when using 6 = 50°. We can clearly 
see that if the majority of the residual emission is lo- 
calised within the central region and is dominated by syn- 
chrotron radiation from an independent, possibly exotic, 
source (e.g. such as dark matter), then we expect that 
the spectral indices describing the synchrotron emission 
in the central and peripheral regions to be significantly 
different. 

We formulate our conclusions based upon results from 
a four-template multi-linear fit to the CMB-subtracted 
WMAP data, involving our free- free, dust and two new 
synchrotron templates, defined by an angle which min- 
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FIG. 7: Central (top) and peripheral (bottom) templates of 
soft synchrotron galactic foreground emission for = 50°, 
presented in units [T(K)], where T is antenna temperature. 




FIG. 8: Xred. verses for 4-template multi- linear fits to the 
CMB-subtracted WMAP data, involving our free- free, dust, 
central and peripheral synchrotron templates when using the 
Gibbs (black) and ILC (magenta) CMB estimators. 



imises the value of Xred (^) associated with the fit. 

In fig. [s] we display the results for x^ed associated with 
4-template fits for different values of the angle 0. We find 
that the minimum value Xred. min. — 6.379 (9.173), occurs 
for 6 48.5° (46.5°) when using the Gibbs (ILC) CMB es- 
timator. Thus, using our 4-template fit we have obtained 
a 16% (17%) reduction in the value of x^d i"elative to 
that obtained with the 3-template fitting procedure dis- 
cussed in § |III[ This can be observed directly in fig. [9) 
where we display maps of the residual emission r and 
maps of the parameter (fj in the K (top), Ka (middle) 
and Q (bottom) frequency bands following the described 
4-template fit when using the Gibbs (left) and ILC (right) 
CMB estimators, where again, for the maps of (fj we 
only display regions of statistically significant emission 
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FIG. 9: Maps of the residual emission (top 3 rows, presented 
in arbitrary units) and of the parameter dj (bottom 3 rows) 
in the K (top), Ka (middle) and Q (bottom) frequency bands 
resulting from our 4-template fit when using the Gibbs (left) 
and ILC (right) CMB estimators, where for the maps of (^fj 



we truncate the maps using a minimum of Q 



(i.e. (fj > 1). The values of the parameter (f for 
the five full-sky residual maps are now 4.43 (5.36), 0.82 
(1.35), 1.01 (2.02), 0.089 (0.17) and 0.028 (0.27) for the 
K, Ka, Q, V and W bands respectively. (The la errors 
on these values are all of order 10~^.) These values cor- 
respond to decreases of 20.00% (18.67%), 7.70% (6.77%), 
6.30% (4.94%), 0.05% (1.40%) and 0.06% (0.35%), rela- 
tive to those obtained for the 3-template fit when using 
the Gibbs (ILC) CMB estimator. (The la errors on these 
values are all of order 10~^ or less.) Again, we compare 
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FIG. 10: Spectral dependences of best fitted soft synchrotron 
galactic foreground emission in the central (green markers) 
and peripheral (magenta markers) regions, resulting from 
our 4-template fit when using the Gibbs (circles) and ILC 
(crosses) CMB estimators. Also plotted are the lines of best 
fit associated with the best fit spectral indices for the dis- 
played data sets (dashed black (Gibbs) and dashed magenta 
(ILC) ). For clarity, we have slightly displaced the frequency 
of the two sets of data markers associated the Gibbs and ILC 
CMB estimators by Alog]^Q(zy/GHz)=H-0.5 and —0.5 respec- 
tively. 



these results for for the full-sky with those calculated 
using a region within 50° of the GC, for which we obtain 
7.94 (9.04), 1.24 (1.81), 1.20 (2.21), 0.107 (0.186) and 
0.027 (0.29) for the K, Ka, Q, V and W bands respec- 
tively when using the Gibbs (ILC) CMB estimator. (The 
la errors on these values are all of order 10~^ or less.) 
These values correspond to decreases in (f of 45.98 ± 
0.12% (45.52 ±0.11%), 24.82 ±0.15% (24.89 ± 0.12%), 
25.24±0.16% (22.00±0.10%), 0.23±0.01% (7.81±0.08%) 
and 0.50 ± 0.01% (2.17 ± 0.01%) in the K, Ka, Q, V and 
W bands respectively, relative to the corresponding 3- 
template fit results. 



In fig. 10 we display the spectral dependences of the 
of the fitted central (green markers) and peripheral (ma- 
genta markers) synchrotron emission components when 
using the Gibbs (circles) and ILC (crosses) CMB es- 
timators. Again, for clarity, we have slightly dis- 
placed the frequency of the two sets of data mark- 
ers associated the Gibbs and ILC CMB estimators by 
Alog]^o(^/G^Hz)=±0.5 and —0.5 respectively. Since the 
differences in the spectral dependence of the two syn- 
chrotron components is very subtle, we have for clar- 
ity omitted plotting curves corresponding to the free-free 
and dust components, which for the current fit, were de- 
duced to possess best fit spectral indices of -1.98 (-1.86) 
for free-free, -2.93 (-2.99) for dust (between K and Ka 
bands) and ±1.14 (±1.23) for (between V and W bands), 
when using the Gibbs (ILC) CMB estimators. 

The respective best fit spectral indices for the cen- 
tral and peripheral synchrotron components were -2.99 
(-3.13) and -2.89 (-3.11). (We have omitted the ILC re- 
sults at 93.5 GHz which gave negative values for the syn- 
chrotron template weights, whose la errors were of the 
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same order, further supporting the decision that such 
results should be deemed unreliable). Both of these val- 
ues are consistent with the range determined in previous 
studies for the spectral index of soft synchrotron emission 
^31^ . Further, both indices are extremely similar in them- 
selves, contrary to the suggestion by Hooper, Finkbeiner 
and Dobler (2007) (HFD henceforth) that the haze emis- 
sion possesses a significantly harder spectral dependence 
(X g] than the conventional soft synchrotron 

emission far from the GC. Also, from fig.[lO] it is clear 
that both regions of synchrotron emission are always 
very similar in intensity, and considering that their spec- 
tral indices are also very similar it is quite possible that 
our results simply indicate the slight regional variation 
in the conventional soft synchrotron emission (i.e. that 
associated with the 408 MHz measurements by Haslam 
et al.)^ rather than evidence for the presence of an inde- 
pendent, exotic source of synchrotron emission near the 
GC. Therefore, we have revealed the interesting result 
that a significant residual haze emission may be repro- 
ducible by simply invoking a slight regional dependence 
in the spectral index of the soft synchrotron foreground 
emission. 



V. SYNCHROTRON RADIATION FROM DARK 
MATTER ANNIHILATION 

One of the most heavily investigated extensions to 
the standard model of particle physics is supersymmetry 
(SUSY), which is a framework based around a proposed 
symmetry between fermions and bosons. Moreover, in 
supersymmetric theories, for every boson (fermion) in the 
standard model there exists a corresponding fermionic 
(bosonic) superpartner that is typically much heavier 
than its standard model counterpart. Such theories have 
great appeal since they can potentially solve the so-called 
hierarchy and unification paradoxes renowned in particle 
physics. (For a brief review of SUSY see [33 §3.2 and 
references therein). 

For our purposes, the greatest appeal of SUSY is that 
it provides an excellent CDM candidate, known as the 
neutralino with properties that lead to a relic density 
^CDM^^ ~ 0.1, consistent with current experimental ev- 
idence [34]. Perhaps most importantly, these properties 
arise independently of the any such astrophysical con- 
straints, but simply from requirements based on their 
ability to solve the hierarchy problem. 

The neutralino is a superposition of higgsinos, winos 
and binos, which are the fermionic superpartners of the 
standard model gauge bosons, and hence results in a di- 
verse range of properties largely unconstrained by present 
observations. Consequently the neutralino is electrically 
neutral and colourless, only interacting weakly and grav- 
itationally, and therefore very difficult to detect directly. 
Hence the neutralino is aptly categorised as a weakly in- 
teracting massive particle or WIMP. In R-parity conserv- 
ing supersymmetric models the lightest neutralino, being 



the lightest SUSY particle (LSP) is stable [35]. Conse- 
quently, in a scenario where present-day CDM exists as 
a result of thermal- freeze out, the dominant species of 
CDM is likely to include the LSP. The relic density of 
the LSP will then heavily depend on its mass and self- 
annihilation cross-section. 

In this section we calculate the expected fiux in 
synchrotron radiation resulting from relativistic elec- 
trons/positrons produced by annihilating neutralino 
dark matter. In order to calculate such fiuxes we need 
several ingredients. Firstly, in § V A| we calculate the en- 
ergy spectrum of injected into the ISM as a result of 
the annihilation of supersymmetric (SUSY) neutralinos 
with described by a specified parameter set. Secondly, 



m 



VB we calculate the steady-state number density of 



, for all relevant energies, at all positions along a par- 
ticular line of site (l.o.s.), when accounting for the effects 
of galactic diffusion. Thirdly, in §|V C| we calculate the 



expected power in synchrotron radiation emitted by each 
electron at the frequency at which we are observing (i.e. 
one of the five WMAP frequencies). Lastly, in 



VD 



we 

use these results to calculate the total integrated power in 
synchrotron radiation along a given l.o.s. for a specified 
model of dark matter. We then compare our estimates 
for the synchrotron fiux for several benchmark neutralino 
dark matter models with our estimates for the residual 
emission resulting from our 3-template fit presented in 



HI 



A. Positron injection spectra of annihilating SUSY 
neutralinos 

We use four benchmark neutralino models that sample 
the vast particle parameter space of the minimal SUSY 
model (MSSM), owing to the fact that neutralinos are 
a superposition of four different supersymmetric parti- 
cles (i.e. winos, binos and two higgsinos). We use four 
benchmark models, the properties of which are as fol- 
lows: Firstly, we consider a 50 GeV neutralino which an- 
nihilates 96% of the time to hh pairs and 4% to r+r~ 
(model 1). Secondly, we consider a 150 GeV gaugino- 
dominated neutralino annihilating with identical branch- 
ing ratios and modes as model 1 (model 2). Thirdly, 
we consider a 600 GeV gaugino-dominated neutralinos 
that annihilate 87% to hh and 13% to r+r~ (model 3). 
Lastly, we consider a 150 GeV higgsino-dominated neu- 
tralinos, that annihilate 58% to VI^+VI^- and 42% to ZZ 
(model 4). 



B. Calculation of steady-state positron distribution 

The spectrum of resulting from the particle cascade 
immediately following the annihilation of two dark mat- 
ter particles is subjected to significant modifications as 
these interact with the ISM. To calculate the spec- 
trum of resulting from diffusion to any given point 
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in the galaxy, we solve for the steady-state solution to 
diffusion-loss equation 



d drie 
dt~de 



K{e,T)V 



dne 
de 



d_ 
de 



6(e,r 



dne 
de 



-Q{e,v), 
(7) 



where dUe/dEe is the number density of per unit en- 
ergy, K{Ee^ r) is the spatial diffusion parameter, b{Ee^ r) 
is the rate of energy loss and Q{Ee^ r) is the source term 
which is proportional to the local rate of dark matter 
annihilation. We solve ^ using exactly the same pro- 
cedure as that described in [36^. Therefore, for brevity, 
we omit a discussion of the full procedure here, instead 
simply describing the parameter sets we utilised in our 
analysis. 

For diffusion and energy loss parameters of the form 



i^(e,r) = i^oe" 6(e, r) 



(8) 



the general form of our solutions for the steady-state 
number density of e^, for a dark matter particle mass 
m^, thermally-averaged annihilation cross-section (mul- 
tiplied by velocity) (cTann.'^) and injection spectrum 
(per annihilation) d^/de, given by 



dne 
de 



1 (P^ 

2 Vm 



{(^ann v)rEe ^ J dc D 

e 

/ drV7(r')exp - ^ ^ > ' 

^ 

oo 



de' 



2rr' 
1)^ 



(9) 



where 



1 



(1-a) 

Ee 

IGeV 



D 

F{e,z) 



[4KoTE {u{e) - i^(e'))]'^' and 
L + {-l)"{2Ln- z) 
D 



-L-h(-l)"(2in-z) 



erf 
—erf 



D 



(10) 



The function /(r), where r is the cylindrical radial coor- 
dinate, is equal to our choice of spherical density profile, 
divided by the scale density ps, averaged over the full 
azimuthal dimension of the diffusion zone (convention- 
ally defined to be the region r > 0, |z| < L containing a 



r. 
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FIG. 11: Positron energy distribution for all positrons con- 
tained within spherical shells of thickness Ar — 200 pc located 
at radii 8 (thick solid), 4 (dashed), 2 (dotted), 1 (dot-dashed) 
and 0.4 (thick solid) kiloparsecs from the Galactic centre, 
calculated using (i) 50 GeV neutralinos with a monochro- 
matic positron spectrum at 50 GeV (black curves) (ii) 50 GeV 
neutralinos with branching ratios B{xx ~^ bb) = 0.96 and 
B{xX ~^ r~^r~) — 0.04 (magenta curves). 



homogeneous and isotropic magnetic field) 



2 , ^2U/2X 



= 2L ^H^^ 



2 , ^2U/2X « 



(r^ + z^) 



-2(/3-7) 



(11) 

As an example, we illustrate the characteristics of the 
steady-state solution by using the simplistic set of param- 
eters similar to that utilised by Finkbeiner (i.e. model 1 
of [2,) ^37, : = i^o = 9 X l{)^^cm^s-^, b = r^^ = 
10~^^s~''^. We use a spherical NFW dark matter density 

with a scale ra- 

25 r-m ~3 



(a,/3,7) = (1,3,1), 



profile [38", 39 , i.e. 

dius and density = 20 kpc and ps = 5.6 x 10~"^ g cm 
respectively. Finally, we use a 50 GeV dark matter par- 
ticle with annihilation cross-section given by ((Jann.'^) = 
2 X 10~2^cm^s~^ which self-annihilates entirely to e+e~. 
This value of ((Jann.^) is very similar to the canonical 
value of 3 X lO^^^cm^s"^ consistent with dark matter 
relic density estimates (^^cdm^^ = 0.1099 ±0.0062) from 
we display the resulting num- 



WMAP 



In fig. 
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ber density as a function of energy, integrated over shells 
of thickness Ar = 200 pc, plotted at radii r =8 (thick 
solid), 4 (dashed), 2 (dotted), 1 (dot-dashed) and 0.4 
(thick solid) kiloparsecs (black curves). These results are 
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clearly consistent with those presented in [2]. In fig. 
we also display the corresponding results for our bench- 
mark model 1, that is, a 50 GeV Gaugino-dominated neu- 
tralino with branching modes B{xX ~^ ^^) = 0.96 and 
B{xX T+r~) = 0.04 and annihilation cross-section 
{o'ann.'^) = 2.7 X lO^^^cm^s"""^ (magcuta curves). 

Comparing the two sets of results we observe that the 
injection spectrum of positrons resulting from dark mat- 
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ter annihilation plays a significant role in the correspond- 
ing steady-state distribution throughout the Galaxy. As 
we expect, the more realistic steeply declining (with re- 
spect to energy), injection spectrum from our model 1 
neutralinos results in a positron distribution which also 
steeply declines with energy at a rate much faster than 
that associated with neutralinos with a monochromatic 
injection spectrum. Such observations prove extremely 
significant when we consider the synchrotron radiation 
fiux produced by these particles since in the following 
subsection we will demonstrate that only above a cer- 
tain energy threshold effectively contribute to the syn- 
chrotron fiux at any given frequency. 



C. Synchrotron radiation spectrum 

Synchrotron radiation from results from their ac- 
celeration by the Galactic magnetic field that defines 
the diffusion zone. For possessing a Lorentz fac- 
tor 7 = E/rrie, confined to helical trajectories with an 
isotropical distribution of pitch angles a, with respect to 
a local magnetic fiux 5, the average fractional power in 
synchrotron radiation emitted at frequency u is given by 
iQlIlT]: 



1 dpsync(^) 27^3 



ptot. 

sync. 



dx 



327r 



da J dyK^iM. (12) 




FIG. 12: The expected synchrotron power to eminate from a 
spherical shell of thickness 200 pc and inner radius 8kpc, at 
the K-band frequency of 22.8 GHz for our realistic neutralino 
model (magenta curve) less-realistic model involving direct 
annihilation channels to electron/positrons (black curve). 



where x = v/ {yBl'^)^ with vb = 3eB / {ATrrriec) :^ 
42 (5/10/iG)Hz, and the total synchrotron power Pgync. 
is set equal to a fraction /sync of the total emitted power 
He) 



a significantly smaller contribution (by at least an or- 
der of magnitude) from higher energy compared to 
the monochromatic produced by our alternative neu- 
tralino model. 



sync. 



/sync.^(e) 
^2 



sync. 



(13) 



where Tgync. 
time scale. 



Te/ f sync. IS the synchrotron energy loss 



In fig. 12 we plot the spectrum in fractional power ( 12 ) 
radiated by electrons in the presence of a magnetic field 
of strength B = 10 /iG at frequency 22.8 GHz (K-band). 
Only electrons with energies above a certain threshold 
7^ > v/vb — 12GeV (i.e. x < 1), contribute efficiently 
to the radiated power. To illustrate the significance of 
this in our present analysis, in fig.[l2] (bottom panel) we 
display the emissivity (i.e. number density multiplied 
by emitted power) in synchrotron radiation at 22.8 GHz 
for a Galactocentric radius of 8kpc, for the monochro- 
matic 50 GeV neutralino dark matter (black line) and the 
50 GeV neutralino dark matter described by our model 1 
(magenta line). Despite the fact that the model 1 neu- 
tralinos generate a proportionally larger number of low 
energy e^, which effectively do not contribute to the syn- 
chrotron power, its steeper injection spectrum results in 



D. Results for synchrotron flux 

We now have all the necessary ingredients in order to 
calculate the synchrotron flux resulting from annihilation 
of Galactic dark matter along a given l.o.s.. To obtain 
such a flux, we flrstly calculate the column density spec- 
trum of e^, da/de by integrating (|9| along a specifled 
l.o.s. 



d(j(J, e) 
de 



dl 



drie 
de 



(14) 



where /max. = L/\sm{S)\ is the path length along our 
l.o.s., which lies within a plane intersecting the GC and 
with inclination S relative to the GC. The spectrum in 
synchrotron power dPgync./di^ from the full column of 
along a given l.o.s. is then simply (in power per area per 
unit solid angle per frequency interval) obtained by inte- 
grating (over energy) the column depth ( 14) with the syn- 
chrotron power spectrum for a single electron, dpsync./di^. 
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given by (12), 



47r 



27^3 



dx du de 

2 dcr((5, e) X 



1287r2 Tsync.^B 



de 



de 



TT CXD 



X X 



cc/sino; 



(15) 



So that we can compare our results with those obtained 
using the treatment by Hooper, Finkbeiner & Dobler 
[8], we use the fohowing set of astrophysical parame- 
ters: Ko = 102^cm3s-\ a = 0.33, te = 2 x lO^^s, 
/sync. = 0.25, B = 10 /iG, L = 3kpc. In fig. 13 we com- 
pare our estimates for the angular distribution of syn- 
chrotron power in the K-band resulting from the annihi- 
lation of dark matter (solid curves), for all four bench- 
mark models, with that obtained from our 3 template- 
based fit to the CMB-subtracted WMAP data when us- 
ing the Gibbs (black markers) and ILC (magenta mark- 
ers) CMB estimators, where once again we have slightly 
displaced the frequency of the two sets of data markers. 
We have also plotted the lines of best fit to the deduced 
K-band residual emission (dashed curves). The noise in 
the measurements of the sky temperature is very small 
at the resolution that we use here. In consequence the 
la statistical error bars on CMB measurements associ- 
ated with the displayed residual emission are also very 
small. Therefore, in order to illustrate these errors we 
have displayed the la error bars, deduced in the proce- 
dure described above, enlarged by a factor of 10. We 
have multiplied each dark matter flux by a universal nor- 
malising factor, or boost factor BF^ in order to obtain 
the best-fit to the WMAP data, where for each model we 
minimise the parameter 



BF x( 

V /DM 



haze- 



a- 



(16) 
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FIG. 13: Estimated synchrotron emission from dark matter 

where {dP^yncXO) / dv)^^^^ and a^(6>) are the differential annihilation using our four benchmark models 1-4 (top to 

residual synchrotron power and (la) error map (squared) bottom respectively), re-normalised, via a boost factor BF, 

in the CMB-subtracted data in the K-band, averaged to the deduced residual emission in the K-band. 
over the beam FWHM at Galacto-centric inclination 9. 

The fiux produced from annihilations within a smooth 
distribution of Galactic dark matter is far less than that 

required to explain the WMAP haze. This is contrary boost factors of approximately 372 (366), 675 (666), 1701 

to the results presented by HFD as well as Caceres & (1680) and 503 (498) for models 1-4 respectively, in order 

Hooper [5 , who claim that such annihilations provide a to generate the necessary synchrotron power and explain 

fiux, which is of the order of several kJy at the relevant our discerned haze emission in the K-band. 

investigated, and therefore fits extremely well to the Despite this, the angular distribution of the estimated 

residual haze emission. Here, we find that we require emission from dark matter fits very well to the data (de- 
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spite the cylindrical nature of our selected dark matter 
density profiles). Interestingly we see that the angular 
dependence of the emission is very similar for all four neu- 
tralino models, with values of equal to 94.4% (98.5%), 
88.7% (92.9%) and 82.1% (86.4%) of that associated with 
model 1, for models 2, 3 and 4 respectively. This implies 
that the gaugino or higgsino fractions are quite irrelevant 
in the determination of the observed synchrotron power 
spectrum. Moreover, the values of the best-fit boost fac- 
tors obtained support the notion that it is the mass of 
the neutralino which dictates the overall normalisation 
of the spectrum (i.e. we find smaller boost factors for 
lighter neutralinos, which possess a higher rate of anni- 
hilation since the annihilation rate scales as ^^^), whilst 
the dark matter profile of the Galaxy dictates the shape 
of the spectrum. 

VI. CONCLUSIONS 

We have investigated the prospects for the existence 
of the WMAP haze and the potential for it to be ex- 
plained with either a regional dependence in soft syn- 
chrotron foreground emission or synchrotron emission re- 
sulting from the diffusion of electrons and positrons pro- 
duced by annihilating neutralino dark matter, and found 
that the former is a better explanation. We conducted 
parallel analyses using both Gibbs and ILC CMB esti- 
mators and found that generally all our results presented 
are mostly insensitive to the choice of the Gibbs or ILC 
CMB estimator, but we believe that the Gibbs estimator 
is somewhat more trustworthy than the ILC. 

Using only a simple 4-template fit, with the syn- 
chrotron template cut into two independent regions with 
different spectral indices, we obtained a significant reduc- 
tion in residual K-band emission compared to our initial 
3-template fit. The reduction was approximately 20% 
(19%) for our full sky analysis and 46% (46%) in the in- 
ner 50° surrounding the galactic centre, when using the 
Gibbs (ILC) CMB estimator. We also found that the 
spectral dependences of the central and peripheral soft 
synchrotron emission possessed very similar spectral de- 
pendences in order to produce this improved fit, implying 
that a small regional variation in the spectral dependence 
of the soft synchrotron foreground emission can account 



for a substantial proportion of the observed residual emis- 
sion. 

Using four benchmark models of SUSY neutralinos we 
calculated the flux in synchrotron radiation in the K- 
band (where the haze possesses the greatest statistical 
significance) resulting from energy losses by electrons and 
positrons produced in dark matter annihilations as they 
diffuse through the ISM. We fitted the angular depen- 
dence of the l.o.s. flux with that of the residual emission 
determined from our 3-template fit. We deduced that 
in all four models, significant enhancements in the an- 
nihilation rate were required in order to account for the 
observed residual emission. Specifically, the calculated 
synchrotron flux resulting from dark matter annihilations 
needed to be multiplied by a boost factor ranging from 
approximately 300-1700 in order to best-fit to the resid- 
ual emission. 

The possibility that such boosting might be achiev- 
able when acknowledging the presence of dark matter 
substructures, which have been postulated to exist to 
various degrees by many recent N-body simulations of 
dark matter halos [42]. This has been the subject of re- 
cent investigation by Borriello, Cuoco and Miele (BCM) 
[43^. BCM conclude that even with such boosting con- 
servative constraints on the annihilation cross-section are 
{o'ann.v) ~ 10~^^ cm^ s~^ — 10~^^ cm^ s~^ , which is up to 
two orders of magnitude larger than the canonical value 
of 3 X 10~^^cm^s~^ consistent with dark matter relic 
density estimates from WMAP. 
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